Strong Noise Effects in one-dimensional Neutral Populations 
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The dynamics of well-mixed biological populations is usually studied by mean-field methods and 
weak-noise expansions. Similar methods have been applied also in spatially extended problems, 
relying on the fact that these populations are organized in colonies with a large local density of 
individuals. We provide a counterexample discussing a one-dimensional neutral population with 
negative frequency-dependent selection. The system exhibits a continuous phase transition between 
genetic fixation and coexistence unexpected from weak-noise arguments. We show that the behavior 
is a non-perturbative effect of the internal noise that is amplified by presence of spatial correlations 
(strong-noise regime). 
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The dynamics of biological populations is traditionally 
modeled using systems of partial differential equations 
for a set of density fields, representing the (local) con- 
centrations of species, resources, genetic traits, etc. pQ. 
Despite the success of deterministic theories, the fate of 
real populations is often determined by purely stochastic 
forces induced by the discreteness of individuals (inter- 
nal noise). In well- mixed populations, the internal noise 
is a finite-size effect that can be studied using system's 
size expansions 2J. Even though the noise has macro- 
scopic effects, including gigantic oscillations and unex- 
pected extinction events [3J, these results are obtained 
in a weak-noise regime accessible from the determinis- 
tic mean-field description by perturbative methods. In 
spatially extended populations, on the contrary, correla- 
tions induced by spatial constraints can enhance the ef- 
fects of the internal noise and induce phenomena that are 
completely unpredictable by means of perturbative meth- 
ods. Spatial populations are usually studied using "patch 
dynamics" [I] or "stepping stone models" |5 , in which 
colonies composed of a number (f2 3> 1) of individuals 
occupy the sites of a lattice and interact by migration. 
All processes inside the colonies (or demes) and between 
them are stochastic, therefore as long as their internal 
population is hnite and mobility is low, the macroscopic 
behavior is affected by noise. Unlike the well-mixed case, 
these phenomena can persist also in the limit of infinitely 
large systems. The simplest situation occurs in neutral 
theories [6], in which purely stochastic genetic (ecologi- 
cal) drift is the only force determining the evolution of the 
whole population, favoring local exclusion and fixation. 
In low dimensional systems, genetic drift competes with 
diffusion and migration processes that instead promote 
genetic mixing [5]. When the microscopic interactions 



are linear, like in the voter model or in the Moran pro- 
cess, changing microscopic details of the model (e.g. il) 
only produces a temporal rescale of the dynamics. In the 
presence of non-linear interactions, instead, the contin- 
uum theory may depend on the microscopic model and 
on the way in which the infinite size limit is taken. 

In this Letter we address this problem considering the 
effects of negative frequency- dependent selection (NFS) in 
a neutral population composed of two competing genetic 
types. In a regime of strong competition, individuals car- 
rying a rare phenotype are able to exploit new resources, 
survive predation or better adapt to an evolving envi- 
ronment, gaining a fitness advantage. This advantage 
is transient, because it declines when the phenotype be- 
comes more common. In this sense, it is a special form of 
Darwinian selection that, within the framework of neu- 
tral theories, favors advantageous mutations and main- 
tains balanced polymorphism and biodiversity. Negative 
frequency-dependent selection has been recently observed 
experimentally in natural populations both at the phe- 
notype and genetic level [7j. While NFS is enough to 
maintain a stable coexistence of different genetic types in 
well-mixed populations, we show that in one-dimensional 
systems the competition between NFS and genetic drift 
leads to an unexpected continuous phase transition be- 
tween genetic fixation and coexistence. 

We consider an idealized population of identical 
haploid individuals that differ in a single gene, i.e. at a 
given locus they have either allele A or B. According to 
Kimura's stepping stone models [S], the population is 
spatially organized in "demes", i.e. islands or colonics 
that occupy the sites of a d-dimensional lattice. Each 
deme contains a fixed number O of individuals. The 
evolution of the population is governed by three stochas- 
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tic processes: 1) a migration process (at rate Ad) that 
consists in the pair-exchange of individuals between 
neighboring denies; 2) a Moran process (at rate Am) in 
which an individual (either of type A or B) is chosen for 
reproduction and her offspring replaces another random 
individual of the same deme; 3) a process of frequency- 
dependent selection at rate As- The latter consists in 
randomly choosing an individual together with other two 
individuals of the same deme. The first individual is re- 
placed by a new one carrying the less/most popular allele 
appearing in the triple. In the case of NFS, less frequent 
alleles are selected. The state of the system is univocally 
described by a vector {n\, n%, . . . , n^a} = n where rii 
is the number of individuals of type A in deme i and 
L is the linear size of the lattice. The evolution of this 
state is governed by the master equation dP(n,t)/dt = 
E„^„ [W(ll' -> n)P(n', t) - W{n -> n')P(n, t)}, with 
local transition rates 

w D (ni -¥ rii - l,rij -> rij + 1) = Ifff— ^ (la) 

v> M (m -> m - l) = A„^ (lb) 

w s (n l ~>n l ~l)= A s (^) 2 ^.(lc) 

Following standard stochastic methods [H |3J, we write 
rii = VLfi + O(Qs), and exploiting the fact that frequen- 
cies change of C(i7 _1 ) in a time step At, we expand 
the master equation in powers of 1/fi. In the limit of 
large fi, after temporal rescaling t = QAt, we obtain a 
Fokker-Planck equation and the corresponding non-linear 
Langevin equation for the local frequency fi of type A. 
We take a naive continuum limit sending the lattice spac- 
ing a to zero and keeping constant the diffusion coefficient 
D = a 2 Ad- The effective continuum equation for f(x,t) 
reads 

| = dvV-K2/-i)/(i-/) + ^/7(i^ (2) 

where v = As, c 2 = (Am + As)/f2 and n(x,t) is a delta 
correlated gaussian noise ({t)(x, t)rj(x', t')) — 8{t—t')5(x— 
x')). The parameter v can be made negative admitting 
also positive frequency-dependent selection. In this case, 
the two symmetric absorbing states at / = 0, 1 are at- 
tractive and the system fixates to one of them after a 
coarsening process. For v > instead, the attractive 
fixed point is the active homogeneous equilibrium state 
at f eq = 1/2, corresponding to a completely mixed co- 
existence. The fluctuations around f eq are gaussian with 
amplitude ©(fi- 1 / 2 ). The mean-field picture, with a dis- 
continuous phase transition at v c = 0, is not correct in 
one dimension. This is evident from the results of nu- 
merical simulations in Fig{T] where we plot the average 
heterozygosity h(t) = L^ 1 (f(x,t)(l — f(x,t)))dx as 
a function of time for a system of size L = 10 3 (with 
Ad = 0.2, Am = 1) and different values of As and fL In a 
finite system, the dynamics cannot escape the absorbing 
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FIG. 1: Time dependence of the average heterozygosity h(t) 
for a system of size L = 10 3 , Ao = 0.2, Am = 1, and different 
values of As- (a) Q. = 5 and A s = 0,0.2,0.22,0.25,0.3. (b) 
= 20 and A s = 0, 0.001, 0.01, 0.02, 0.05. Numerical data are 
averaged over 100 initial conditions. 

state and h eventually vanishes. However, the temporal 
behavior of h clearly shows the existence of a phase tran- 
sition between fixation and coexistence at some non-zero 
As (i.e. v c > 0). For il = 5 the coexistence phase de- 
velops gradually, indicating the continuous nature of the 
transition. Increasing fl (Figjlja), the critical point moves 
closer to zero, and the transition becomes more abrupt. 
In this case, numerical results do not allow to distinguish 
between continuous or discontinuous behavior. A discon- 
tinuous transition is recovered for sure in the thermody- 
namic limit, at least if we first take Q — > oo and then 
L — > oo limits. An argument in favor of the continuous 
character of the transition is based on Hinrichsen's con- 
jecture about the impossibility of first-order transitions in 
non-equilibrium Id systems with short-range interactions 
(5]. Indeed, in the absence of surface tension (d — 1), a 
small droplet of the absorbing phase can always grow 
by diffusion of the boundaries and eventually destabilize 
the fluctuating active equilibrium. For fl < oo, the nucle- 
ation rate Re of a droplet of size I of the absorbing phase 
out of the homogeneous equilibrium is small but finite 
(Re ~ e - c ™ with c oc ^ +\ )• Hence, if the only pro- 
cesses into play were nucleation and diffusion, the absorb- 
ing phase would always dominate. It is thus necessary to 
have a different mechanism for the emergence of coex- 
istence. According to this argument, the Fokker-Planck 
approach and small-noise expansions are not valid even in 
the large il limit, calling for a different "non-perturbative 
approach" . Interestingly, a non-perturbative renormal- 
ization group study of Eq[2]was done in [9] showing that, 
for d < d' c fa 4/3, the critical behavior of the model is 
controlled by a non-Gaussian fixed point of the RG flow. 
In fact, [5] suggests the existence of a continuous phase 
transition at v c > that cannot be found by expanding 
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in I/O, around the mean-field theory. 

A way to confirm and understand this picture is that of 
studying a model with only one individual per site and in- 
teractions with neighboring sites. The migration process 
is unchanged and consists of exchanging two neighboring 
individuals with rate Xd- The Moran process is reduced 
to a voter-like move between nearest neighbors. With 
rate Am, an individual dies and is immediately replaced 
by an offspring of a randomly selected neighbor. We 
adopt a spatial version of NFS recently proposed in the- 
oretical ecology to model interspecies competition [TU] , 
Two neighboring individuals that are different compete 
for reproduction: with a rate As we randomly select a 



neighbor of the two to be replaced with the less com- 
mon type among the three selected individuals. In the 
absence of diffusion (Ad = 0), the one-dimensional dy- 
namics can be mapped on that of branching-annihilation 
random walks with even offsprings (BARWe) [TT] or it 
can be analyzed with approximated techniques developed 
for binary spin systems. After the mapping <7i = 2ri^ — 1, 
the transition rates take the simple form WM{<Ji ~> 

the master equation for P(<j, £), we can compute coupled 
evolution equations for the multi-spin correlation func- 
tions {(Ji . . . tjj) that cannot be solved exactly. At the 
lowest order of the hierarchy we have 



= ~2~ Wi-V + W+l) ~ 2 \ a i)\ + ~2 WiVi-^i-l) + (< J i< J i+l (J i+2) ~ 2(0i)] (3) 



Some progress can be done by assuming the Kirkwood 
approximation, i.e. factorizing three-sites and four-sites 
quantites as products of lower order ones [12] . For con- 
venience we define the average frequency of type A as 
/, = (1 + (cr,-))/2 and the average local heterozygosity 
hi = (1 — (o"j(Ti+i))/2. Taking the continuum limit, and 
neglecting higher order contributions, we get equations 
for f(x,t) and h(x,t), 

%=°S-Xs(2f-l)h (5) 
dh d 2 h 

°g - (D M +D s )^-2h(\ M -\ s (l-2h)) (6) 

where Dm = cl 2 Xm and Ds(h) = a 2 \sh{x,t). The ef- 
fective diffusion terms generated in the continuum limit 
account for the random walk of domain walls induced 
by both Moran and NFS moves. The coefficient Dg(h) 
depends linearly on the local heterozygosity h(x, t) as 
required by NFS microscopic rule. Eqs. [5][6] are deter- 
ministic and describe the average behavior of the system. 
A single realization of the process would also present a 
noise term with two homogeneous absorbing states at 
f(x,t) = 0, 1 Vx. The stationary state is obtained solv- 
ing Eqs. [5][6]with the assumption of spatial translation 
invariance, f(x) ~ f and h(x) ~ h. There are two sta- 
tionary solutions of Eq{5] depending on the ratio As /A m- 
A continuous phase transition separates the absorbing 
phase (h — 0) where a unique species survives from an 
active phase of non-trivial coexistence (h = 4^ with 
A = As/Aa/). The theoretical curve Ji(Xs/Xm) is plotted 
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in Fig{2] (full line) together with the results of numeri- 
cal simulations (points). For Xd = the agreement is 
very good, although in the analytical calculation we ne- 
glected the effects of the self-generated spatial diffusion 
that could be responsible of the renormalization of both 
the critical point and the critical exponents. To explore 
the dynamics in the absorbing phase, we prepared the 
population of size L = 10 in a completely mixed con- 
figuration (/ = 0, h = 1/2) and we studied the tem- 
poral behavior of the average cluster size (of type A or 




FIG. 2: Average heterozygosity h in the stationary state as a 
function of the control parameter Ag (with Am = 1) for the 
model with single individual per site. Numerical simulations 
(points) are performed on a system of size L — 10 3 for differ- 
ent migration rates: Xd — (blu circles), 2 (red diamonds), 
20 (green squares). The case A_d = is in good agreement 
with the stationary solution of Eqs. |5)6| (full line). 
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FIG. 3: (a) 1 + 1 dimensional snapshot of the evolution of a 
system of size L = 10 4 starting from random initial conditions 
into the absorbing phase (As = 0.5); (b) temporal decay of the 
average local heterozygosity h for As = 0.5, the dashed line 
indicates a power-law decay oc t -0,5 ; (c-d) 1+1 dimensional 
snapshots of the evolution of a system of size L = 10 4 starting 
from a single defect towards the stationary active phase for 
A s = 1.05 (c) and 1.5 (d). 

B) until fixation. Clusters growth is evident from the 
time-oriented snapshots in Fig. [3^,. The average domain 
length is roughly lit) oc l/h(t), therefore the behavior 
Eq. [6] provides information on the coarsening exponent. 
In the absorbing region, we can neglect the non-linear 
terms in the r.h.s. and solving the linear equation we ob- 
tain a power law decay h(t) oc i -1 / 2 . This is confirmed 
by the numerical results shown in Fig. [3]d. The typi- 
cal domain length £(t) thus grows in time as £(t) <~ t 1 / 2 , 
which is the same behavior as the pure voter model (i.e. 
As = 0). This law is not modified by small diffusion, 
but the case VL > 1 requires further investigation (see 
FigjTJ) . In the active phase, the system shows local clus- 
tering up to a characteristic length scale £ oc 1/h, that 
is reduced by increasing As from just above criticality 
(Fig{3j;) to the bulk of the coexistence region (FigpH). 
Comparing snapshots in Figs|3]>[3ji with that in Fig]3^,, 
we notice a much higher local volatility of patches and 
domain walls, a signature of increasing genetic mixing. 
The addition of a small migration rate does not change 
the qualitative picture observed for Xjj = 0. However, in- 
creasing diffusion favors local mixing of individuals and 
acts against genetic drift. This is confirmed by numeri- 
cal simulations reported in Fig[2] that show the critical 
point moving towards zero and the curve of h becoming 
steeper for increasing values of A^i. In the limit A^> — > oo 
the constraints imposed by the one-dimensional geometry 
arc eliminated leaving a completely well-mixed popula- 
tion that should undergo a first-order transition at zero 
NFS as predicted by the mean-field theory. 

In conclusion, negative frequency-dependent selection 
is a major force for coexistence, but spatial constraints 



tend to amplify the effects of genetic drift promoting fix- 
ation. The outcome of this competition is a phase tran- 
sition between a monomorphic and a polymorphic phase. 
In one spatial dimension, the nature of the phase tran- 
sition critically depends on the number Q of individuals 
per lattice site and on their mobility Ad- In the infi- 
nite size limit, any finite value of f2 generates a dynamics 
that cannot be obtained from the mean-field by means 
of perturbative approaches (strong noise regime). Re- 
cent works have shown that the strong noise limit can be 
important in the case of weak directional selection, there- 
fore the present results could be applied to a large class 
of phenomena occurring at the one-dimensional fronts of 
expanding populations, such as Fisher's waves [13] and 
"gene surfing" [2]. Negative frequency-dependent selec- 
tion is also expected to play a role in plant communities 
and marine ecosystems. It would be interesting to inves- 
tigate whether one-dimensional structures such as val- 
leys, riverbanks and seacoasts, could act as spatial bot- 
tlenecks promoting or impeding genetic variability even 
in higher dimensional populations. 
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